Randomized Gram--Schmidt Process with Application to GMRES
نویسندگان
چکیده
A randomized Gram--Schmidt algorithm is developed for orthonormalization of high-dimensional vectors or QR factorization. The proposed process can be less computationally expensive than the classical while being at least as numerically stable modified process. Our approach based on random sketching, which a dimension reduction technique consisting in estimation inner products by their small efficiently computable images, so-called sketches. In this way, an approximate orthogonality full obtained orthogonalization provide computational cost any architecture. benefit sketching amplified performing nondominant operations higher precision. case numerical stability guaranteed with working unit roundoff independent problem. applied to Arnoldi iteration and results new Krylov subspace methods solving systems equations eigenvalue problems. Among them we chose GMRES method practical application methodology.
منابع مشابه
Gram–Schmidt Process in Different Parallel Platforms
Important operations in numerical computing are vector orthogonalization. One of the well-known algorithms for vector orthogonalisation is Gram–Schmidt algorithm. This is a method for constructing a set of orthogonal vectors in an inner product space, most commonly the Euclidean space Rn. This process takes a finite, linearly independent set S = {b1, b2, ..., bk} vectors for k ≤ n and generates...
متن کاملA Parallel Variant of the Gram-Schmidt Process with Reorthogonalization
c © 2006 by John von Neumann Institute for Computing Permission to make digital or hard copies of portions of this work for personal or classroom use is granted provided that the copies are not made or distributed for profit or commercial advantage and that copies bear this notice and the full citation on the first page. To copy otherwise requires prior specific permission by the publisher ment...
متن کاملModified Gram-Schmidt (MGS), Least Squares, and Backward Stability of MGS-GMRES
The generalized minimum residual method (GMRES) [Y. Saad and M. Schultz, SIAM J. Sci. Statist. Comput., 7 (1986), pp. 856–869] for solving linear systems Ax = b is implemented as a sequence of least squares problems involving Krylov subspaces of increasing dimensions. The most usual implementation is Modified Gram-Schmidt GMRES (MGS-GMRES). Here we show that MGS-GMRES is backward stable. The re...
متن کاملRounding error analysis of the classical Gram-Schmidt orthogonalization process
1 CERFACS, 42 Avenue Gaspard Coriolis, 31057 Toulouse Cedex 1, France ([email protected]). 2 The University of Tennessee, Department of Computer Science, 1122 Volunteer Blvd., Knoxville, TN 37996-3450, USA ([email protected]). 3 Institute of Computer Science, Academy of Sciences of the Czech Republic, Pod vodárenskou věž́ı 2, CZ-182 07 Prague 8, Czech Republic ([email protected]). 4 Heinrich-Heine-...
متن کاملSVP , Gram - Schmidt , LLL
Last time we defined the minimum distance λ1(L) of a lattice L, and showed that it is upper bounded by √ n · det(L)1/n (Minkowski’s theorem), but this bound is often very loose. Some natural computational questions are: given a lattice (specified by some arbitrary basis), can we compute its minimum distance? Can we find a vector that achieves this distance? Can we find good approximations to th...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
ژورنال
عنوان ژورنال: SIAM Journal on Scientific Computing
سال: 2022
ISSN: ['1095-7197', '1064-8275']
DOI: https://doi.org/10.1137/20m138870x